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We develop a dynamic point process model of correlated default 
timing in a portfolio of firms, and analyze typical default profiles in 
the limit as the size of the pool grows. In our model, a firm defaults 
at a stochastic intensity that is influenced by an idiosyncratic risk 
process, a systematic risk process common to all firms, and past de- 
faults. We prove a law of large numbers for the default rate in the 
pool, which describes the "typical" behavior of defaults. 

1. Introduction. The financial crisis of 2007-09 has made clear the need 
to better understand the diversification of risk in financial systems with 
interacting entities. Prior to the crisis, the common belief was that risk had 
been diversified away by using the tools of structured finance. As it turned 
out, the correlation between assets was much larger than supposed. The 
collapse fed on itself and created a spiral. 

We study the behavior of defaults in a large portfolio of interacting firms. 
We develop a dynamic point process model of correlated default timing, 
and then analyze typical default profiles in the limit as the number of con- 
stituent firms grows. Our empirically motivated model incorporates two dis- 
tinct sources of default clustering. First, the firms are exposed to a risk factor 
process that is common to all entities in the pool. Variations in this system- 
atic risk factor generate correlated movements in firms' conditional default 
probabilities. Das, Duffie, Kapadia and Saita [5] show that this mechanism 
is responsible for a large amount of corporate default clustering in the U.S. 
Second, a default has a contagious impact on the health of other firms. This 
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impact fades away with time. Azizpour, Giesecke and Schwenkler [1] pro- 
vide statistical evidence for the presence of such self-exciting effects in U.S. 
corporate defaults, after controlling for the exposure of firms to systematic 
risk factors. 

More precisely, we assume that a firm defaults at an intensity, or con- 
ditional arrival rate, which follows a mean-reverting jump-diffusion process 
that is driven by several terms. The first term, a square root diffusion, rep- 
resents an independent, firm-specific source of risk. The second term is a 
systematic risk factor that influences all firms, and that generates diffusive 
correlation between the intensities. For simplicity, we take this systematic 
risk factor to be an Ornstein-Uhlenbeck process. The third term affecting 
the intensity is the default rate in the pool. Defaults cause jumps in the in- 
tensity; they are common to all surviving firms. We thus have two sources of 
correlation between the firms: the dependence on the systematic risk factor 
and the influence of past defaults. While this formulation parsimoniously 
captures several of the sources of default correlation identified in empirical 
research, the intricate event dependence structure presents a challenge for 
the mathematical analysis of the system. 

Our goal is to understand the behavior of the default rate in the portfo- 
lio in the limit as the number of firms in the pool grows. Large stochastic 
systems often tend to have macroscopic organization due to limit theorems 
such as the law of large numbers. This allows us to identify typical behavior. 
Our main result is a law of large numbers for the default rate in the pool; 
this describes the macroscopically typical profile. The limiting default rate 
satisfies an integral equation that makes explicit the role of the contagion ex- 
posure for the behavior of default clustering in the pool. The result depends 
heavily on the analysis of Markov processes via the martingale problem; see 
Ethier and Kurtz [11]. We will have more to say on the mathematical as- 
pects of this in a moment. Once the typical behavior has been identified, 
one can then study Gaussian fluctuations and the structure of atypically 
large default clusters in the portfolio. We plan to pursue these directions in 
a future work. 

Previous studies have analyzed the behavior of defaults in large portfo- 
lios. Dembo, Deuschel and Duffie [9] examine a doubly-stochastic model of 
default timing. In their model, default correlation is due to the exposure of 
firms to a common systematic risk factor which is represented by a random 
variable. Conditional on this variable, defaults are independent. A large de- 
viation argument leads to an approximation of the tail of the conditional 
portfolio loss distribution. Glasserman, Kang and Shahabuddin [14] study 
a copula model of default timing using large deviation techniques. In that 
formulation, default events are conditionally independent given a set of com- 
mon risk factors. Bush, Hambly, Haworth, Jin and Reisinger [2] prove a law 
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of large numbers for a related dynamic model. Davis and Rodriguez [6] de- 
velop a law of large numbers and a central limit theorem for the default 
rate in a stochastic network setting, in which firms default independently of 
one another conditional on the realization of a systematic factor governed 
by a finite state Markov chain. Sircar and Zariphopoulou [20] examine large 
portfolio asymptotics for utility indifference valuation of securities exposed 
to the losses in the pool. As with these papers, our model includes exposure 
to a common systematic risk factor. In contrast, however, our model captures 
the self-exciting nature of defaults. Therefore, the firms in the pool are cor- 
related even after conditioning on the path of the systematic factor process. 

The use of interacting particle systems to study the behavior of default 
clustering in large portfolios is a growing area. In a mean-field model, Dai 
Pra, Runggaldier, Sartori and Tolotti [3] and Dai Pra and Tolotti [4] take 
the intensity of a constituent firm as a deterministic function of the percent- 
age portfolio loss due to defaults. In a model with local interaction, Giesecke 
and Weber [13] take the intensity of a constituent firm as a deterministic 
function of the state of the firms in a specified neighborhood of that firm. 
The interacting particle perspective leads to the study of the convergence 
of interacting Markov processes, laws of large numbers for the percentage 
portfolio loss, and Gaussian approximations to the portfolio loss distribution 
based on central limit theorems. The interacting particle system which we 
propose and study incorporates an additional source of clustering, namely, 
the exposure of a firm to a systematic risk factor process. Moreover, firm- 
specific sources of default risk are present in our system. Also, the nature of 
mean-field interaction in our system is different. In [3] and [4], a constituent 
intensity is a function of the current default rate in the pool. In that for- 
mulation, the impact of a default on the dynamics of the surviving firms is 
permanent. In our work, a constituent intensity depends on the path of the 
default rate. The impact of a default on the surviving firms is transient, and 
fades away exponentially with time. There is a recovery effect. 

As we were finishing this work, we learned of a related law of large numbers 
type result by Cvitanic, Ma and Zhang [16]. They take the intensity of a 
constituent firm as a function of a firm-specific risk factor, a systematic risk 
factor and the percentage portfolio loss due to defaults. The risk factors 
follow diffusion processes whose coefficients may depend on the portfolio 
loss. Our model of the risk factors is more specific than theirs, and thus we 
are able to arrive at slightly more explicit results. Moreover, the effect of 
defaults in [16] is permanent, as in [3] and [4]. 

There are several mathematical contributions in our efforts. Our analysis 
of typical events (a weak convergence result) is somewhat similar to that of 
certain genetic models (most notably the Fleming- Viot process; see Chap- 
ter 10 of [11], Fleming and Viot [12] and Dawson and Hochberg [7]), but the 
specific form of our intensity processes imply both complications and simpli- 
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fications. Our work is centered on a jump-diffusion intensity process which is 
driven by Ornstein-Uhlenbeck and square root diffusion terms. This formu- 
iation ailows some explicit simplifying calculations which are not available in 
a more abstract framework. On the other hand, due to the square root sin- 
gularity, certain technical estimates need to be developed from scratch (see 
Section 10). A final point of interest is heterogeneity. Interacting particle 
systems are often assumed to have homogeneous dynamics, where various 
parameters are the same for each particle. This allows the main mathe- 
matical arguments to take their simplest form. Practitioners in credit risk, 
however, in reality face an extra problem in data aggregation, where each 
firm in a portfolio has its own statistical parameters. We have framed our 
weak convergence result to allow for a distribution of "types," that is, a fre- 
quency count of the different model parameters. This leads us to the correct 
effective dynamics of the portfolio and, in particular, to a precise formulation 
of the effects of self-excitation (see Remark 5.2). 

The rest of this paper is organized as follows. Section 2 formulates our 
model of default timing. We establish that our model is well-posed via the 
results of Section 3. In Section 4 we identify the limit as the number of firms 
in the portfolio goes to infinity — a law of large numbers result. The proof 
of this result is in Section 8, but depends upon the technical calculations of 
Sections 5, 6 and 7. Section 9 concludes and discusses extensions. Section 10 
contains a number of technical results on square-root-like processes which 
are used in our calculations. 

2. Model, assumptions and notation. We construct a point process model 
of correlated default timing in a portfolio of firms. We assume that (fi, P) 
is an underlying probability triple on which all random variables are defined. 
Let {W n } n gN be a countable collection of standard Brownian motions. Let 
{ e n}neN be an i.i.d. collection of standard exponential random variables. 
Finally, let V be a standard Brownian motion which is independent of the 
W n, s and e n 's. Each W n will represent a source of risk which is idiosyn- 
cratic to a specific firm. Each t n will represent a normalized default time for 
a specific firm. The process V will drive a systematic risk factor process to 
which all firms are exposed. 

Fix an iV E N, n £ {1,2, ... , N} and consider the following system: 



N,n 




N,n 



\N,n)dt + y X t ,n dW\ 



dL t + ejv/8jv )n A t ' dX t 



t>0, 



(2.1) 



A, 



N,n 





dX t 



-fXtdt + dVt, 



t>0, 
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X = x , 

^ = ^Ex [ e„ ) oo)(£ o A^d S ). 

Here, /3^ n G M + = f [0,oo) and /3}y n S M are constants which represent the 
exposure of the nth firm in the pool to L N and X, respectively. The ckat^'s, 
A/v,n's and (JN,nS are in M + and characterize the dynamics of the firms. 
We will address the role of £n in a moment. The initial condition x of 
X is fixed and 7 > 0. We use x to represent the indicator function here 
and throughout the paper. The description of L N is equivalent to a more 
standard construction. In particular, define 

(2.2) T^ d =inf{i>0:^ Af' n (fe>e n j. 

Then 

( 2 -3) X[t n ,oo) (J Af ds^j = X{ r iv,n< i} 

and, consequently, 

1 N 

n=l 

The process \ N,n represents the intensity, or conditional event rate, of the 
nth firm in a portfolio of N firms. More precisely, X N,n is the instantaneous 
Doob-Meyer compensator to the default indicator process (2.3); see (4.1). 
We will see in Proposition 3.3 in Section 3 that the A^'^'s are indeed non- 
negative. The process X represents a source of systematic risk; in our model 
this is a stable Ornstein-Uhlenbeck process. The process L N is the default 
rate in the pool. The jump-diffusion model for X N,n captures several sources 
of default clustering. A firm's intensity is driven by an idiosyncratic source 
of risk represented by a Brownian motion W n , and a source of systematic 
risk common to all firms — the process X. Movements in X cause correlated 
changes in firms' intensities and thus provide a source of default clustering 
emphasized by [5] for corporate defaults in the U.S. The sensitivity of X N,n 
to changes in X is measured by the parameter fifj n . The second source 
of default clustering is through the feedback ( "contagion" ) term n dL^ . 
A default causes an upward jump of size j}P% n in the intensity \ N > n . Due 

to the mean-reversion of A ,n , the impact of a default fades away with time, 
exponentially with rate ctN tn . Self-exciting effects of this type have been 
found to be an important source of the clustering of defaults in the U.S., 
over and above any clustering caused by the exposure of firms to systematic 
risk factors [1]. 
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In the special case that fi% n = fi% n = for all n e {1, 2, . . . , N}, the inten- 
sities A ,n follow independent square root processes so firms default inde- 
pendently of one another. The formulation (2.1) is a natural generalization 
of the widely used square root model to address the clustering between de- 
faults. 

The interest in large pools of assets is that they provide diversification; 
they allow one to construct portfolios which have small variance. The dy- 
namics of X imply that X is stochastically of order 1, that is, it is stable. 1 
Thus, the only way for the pool to have small variance in our model is for 
each of the constituent firms to have small exposure to X . We thus assume 
that 

lim en = 0. 

If En is not small, the influence of the systematic risk factor X will be of 
order 1, and the "typical" behavior of the pool will strongly depend on X 
(and the tail behavior of the whole system will be strongly determined by 
the tail of X). 

Remark 2.1. Given the simple structure of X, our model is equivalent, 
if x = 0, to a model where each intensity has exposure of order 1 to a small 
systematic risk. Namely, if x = 0, then EjyX = X N where 

dX? = -7-Xf dt + e N dV t . 

Our model allows for a significant amount of bottom-up heterogeneity; the 
intensity dynamics of each firm can be different. We capture these different 
dynamics by defining the "types" 

(2-4) p N ' n = f (aN,n,^N,n,VN,n, Pn >w PN,n)i 

the p^' n 's take values in parameter space V = f x M. In order to expect 
regular macroscopic behavior of as X — y oo, the p^ ,n, s and the A 0; _/v,n's 
should have enough regularity as N — > oo. For each Af£N, define 

1 N 1 N 

7T n = -^2 6 p n -" and A o=jjY, 6 ^^ 

n=X n=l 

these are elements of ^(V) and ^(R-j.), respectively. 2 

We need two main assumptions. First, we assume that the types of (2.4) 
and the initial distributions (the A ,Ar in 's) are sufficiently regular. 



1 Regulatory agencies, for example, are charged with preventing systematic factors from 
spiraling out of control. 

2 As usual, if E is a topological space, &(E) is the collection of Borel probability 
measures on E. 
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Assumption 2.2. We assume that 

7r *== lim ir N and A *== lim A^ 

N-^oo >oo 

exist [in 3 s (V) and ^(M+), resp.]. 

Note that this is what happens in practice; one constructs a frequency 
count of the parameters of the different assets in a large pool and uses this 
to seek aggregate dynamics for the pool itself. For a large pool, one hopes 
that this frequency count will have some simpler macroscopic description. 
Second, we assume that the types are bounded. 

Assumption 2.3. We assume that there is a K 2 .3 > such that the 
"iV.n's, A/v, n 's, CTjv.n's, l$v,J' s > l$v,J' s and ^o,7V,n's are all bounded by K 2 .z 
for all N G N and n G {1, 2, . . . , N}. 

Equivalently, we require that the 7T/v's and A^'s all (uniformly in N) have 
compact support. We could relax this requirement, at the cost of a much 
more careful error analysis. 

We are interested in the typical behavior of {L N }. In Section 3 we consider 
the well-posedness of the model (2.1), while in Section 4 we state the law of 
large numbers result, Theorem 4.2. 

3. Well-posedness of the model. We here state several technical results 
concerning the intensities which are a central part of our model. We want 
to understand the structure of the X ,n, s a bit more. The complications 
which require our attention are the square root singularity, and the fact 
that the XtdXt term contains the term XtXtdt, implying that the dynamics 
of the M 2 -valued process (A,X) contain a superlinear drift. The proofs of 
the results here will be given in Section 10. 

Let W* be a reference Brownian motion with respect to a filtration 
{@t}t>o- Assume also that V is adapted to {^t}t>o- Let £ be a {&t}t>o- 
adapted, point process which takes values in [0, 1] and such that £o = 0. Fix 
p = (a, A, a, /3 C , /3 5 ) G V and A Q in R+. Consider the SDE 



d\ t = -a(X t -X)dt + cr^X t V0dW t * + (3 C d£ t + /3 S X t dX t , t > 0, 

(3.1) 

Ao — A . 

Note that by expanding the dynamics of dX and rearranging a bit, we get 
that 



dX t = -{a + P s -yX t }X t dt + aXdt + /3 C d£ t + a^X t V0dW t * + (3 s X t dV t . 

Also, we have for the moment subsumed the small parameter en into the 
f5 s term; see the proof of Proposition 3.3. 

We will use a number of ideas from [15] (see also [8]). 
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Lemma 3.1. There is a nonnegative solution A of the M-valued SDE 
(3.1). Furthermore, sup t6 [ 0j T] ^[|^*l p ] < 00 f or all T > and p > 1. 

We also have uniqueness. 

Lemma 3.2. The solution of (3.1) is unique. 

The model (2.1) is thus well posed. 

Proposition 3.3. The system (2.1) has a unique solution such that 
Af ' n > for every N G N, n G {1, 2, . . . , N} andt>0. 

Proof. Using Lemmas 3.1 and 3.2, solve (2.1) between the default 
times. Replace (3 s by en(3 s in applying Lemma 3.1. □ 

We shall also need a macroscopic bound on the intensities. 

Lemma 3.4. For each p>l and T > 0, 

1 N 

iW/=' sup ^£ E 0Af' T ] 

is finite. 

4. Typical events: A law of large numbers. Our first task is to under- 
stand the "typical" behavior of our system. To do so, we need to understand 
a system which contains a bit more information than the default rate L . 
For each JV£M and n £ {1, 2, ... , N}, define 

(4.1) Mf'"^ f X[ o,e n )(£ o Af in ^) =X {T ^ >t} 

[where r N,n is as in (2.2)]. In other words, M^' n = 1 if and only if the nth 
firm is still alive at time t; otherwise M^'™ = 0. Thus, M^'" is nonincreasing 
and right-continuous. It is easy to see that 

Mf'"+ f \^ n M^ n ds 
Js=0 

is a mar tingale. Define V = VxR+. For each N € N, define pf = (p^, Af ' n ) 
for all n £ {1, 2, . . . , N} and t>\. For each t > 0, define 
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in other words, we keep track of the empirical distribution of the type and 
credit spread for those assets which are still "alive." We note that 

for all t > 0. 

We want to understand the dynamics of /x^ for large iV (this will then 
imply the "typical" behavior for Lf). To understand what our main result 
is, let's first set up a topological framework to understand convergence of 
fj, . Let E be the collection of sub-probability measures (i.e., defective prob- 
ability measures) on V, that is, E consists of those Borel measures v on V 
such that v(P) < 1. We can topologize E in the usual way (by projecting 
onto the one-point compactification of V\ see [19], Chapter 9.5). In partic- 
ular, fix a point * that is not in V and define V + = f V U {*}. Give V + the 
standard topology; open sets are those which are open subsets of V (with 
its original topology) or complements in V + of closed subsets of V (again, 
in the original topology of V). Define a bijection l from E to £P{P + ) (the 
collection of Borel probability measures on P + ) by setting 

(tu)(A) d ^ f u(A nV) + (l- y{V))8M) 

for all A G 33(V + ). We can define the Skorohod topology on &(P + ), and 
define a corresponding metric on E by requiring i to be an isometry. This 
makes E a Polish space. We thus have that fj, is an element 3 of De[0, oo). 

The main theorem of this section is Theorem 4.2, essentially a law of 
large numbers. The construction of the limiting process will take several 
steps. First, for each p = (a, A, a, (3 C , (3 s ) G V, let b p satisfy 

b p (t) = l-W(b p (t)f -ab p (t), t>0, 

(4.2) 

6 p (0) = 0. 

Note that if b p {t) = 0, then b p (t) = 1 > 0. Thus, b p (t) > for all t > 0. 

The next lemma is essential for the characterization of the limit. Its proof 
is deferred to Section 10. 

Lemma 4.1. There is a unique M + -valued trajectory {Q(t);t > 0} which 
satisfies the equation 

Q{t)= [ P c \b p (t)\+ I b p {t-r){Q{r) + a\}alr 

Jp=_(p,\)eV L Jr=0 

P =(a,\,<r,P°,P S ) 



3 If S is a Polish space, then Ds[0,oo) is the collection of maps from [0, oo) into S 
which are right-continuous and which have left-hand limits. The space Ds[0, oo) can be 
topologized by the Skorohod metric, which we will denote by ds; see Chapter 3.5 of [11]. 
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x exp 



-b p (t)X 



b p {t-r){Q{r) + a\}dr 



r=0 



x 7r(dp)A Q (dA). 
Here, ir and A Q are as in Assumption 2.2. 

Now let W* be a reference Brownian motion. For each p = (p, A ) €. V 
where p = (a, A, a, /3 C , /3 s ), let A£(p) be the unique solution to 



(4.4) 



X* t (p) = X -a[ (X* s (p)-X)ds + a [ ^K¥)dW* s 

Js=0 Js=0 



+ 



Q(s) ds. 



s=0 



We now have our main result. 

Theorem 4.2. For all A £ B§{V) and B e ^(K+), define 

ct 



ik(A x B) d M [ XA (p)E 



(4.5) 

Then 
(4.6) 

for every 5 > 0. Define 
(4.7) 



Xs(A?(p))exp 



A*(p) ds 



s=0 



x 7r(dp)A (dA). 
Km P{^ (73) (/i 7V ,/i)>5} = 



F{t) d Ml- h {V) 



1 



E 



'p=( P ,A)eP 
Then, for all 5 > and T > 0, 



exp 



A*(p)ds 



s=0 



7r(dp)A (dA). 



lim P{ sup \L? -F(t)\>6\=0. 



-0<t<T 

The parts of the proof of this result will be given in Sections 5, 6 and 7. 
In particular, in Section 5 we identify a candidate limit for {// } using the 
martingale problem formulation. Then in Section 6 we prove that {/J, N } is 
tight, which ensures that the laws of {/j, N ys have at least one limit point. In 
Section 7 we prove that the limit is necessarily unique. Then, in Section 8 
we collect things together to prove Theorem 4.2. 

With this result in hand, we can rewrite (4.4) to see the exposure of a 
typical firm to the contagion factor. 
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Remark 4.3. We have that 



F(t) 



(4.8) 



Thus, 



(4.9) 



where 



E 



p=( p ,x)ev 



A f *(p)exp 



A*(p) ds 



s=0 



ir(dp)A (dX) 



P =(p,x)eV 



At(p) = A -a 



f (A* (p) -X)ds + a f dW* 

Js=0 JO 



+ 



B(p s )F(s)ds, 



dcf 



/3 (7 A/i(dp) 



'p=(p,A)G^ 

p=(a,A,<T,^ c ', i 8 s ) 



'p=( P ,A)eP 

P=(q,A,o-,/3 c ',/3 s ) 



A/i(dp) 



for all /j, £ E. In other words, the effective sensitivity of a typical inten- 
sity to the contagion is given by an average weighted by the instantaneous 
intensities. Note that < -B(/x) < i^2.3- 

The homogeneous case provides more explicit insights into the role of the 
contagion exposure for the behavior of default clustering in the pool. 

Remark 4.4. Fix p= (p, A ) G V where p= (a, A, a, f3 c , /3 s ). Assume 
that the pool is homogeneous, that is, p N > n = p for all JVgN and n £ 
{1,2,...,N}. By the relation (4.9), we then have that Q(t) = (3 c F(t). In 
this case, F is given by the unique solution to the integral equation 



F(t) = 1 - exp 



aX / b p (t-r)dr-P 



c 



r=0 



F(r)b p (t-r)dr -b p (t)X c 



r=0 



Furthermore, if there is no exposure to contagion, that is, f3 c = 0, then this 
integral equation reduces to the well-known explicit formula 



F(t) = 1 - exp 



-aX 



b p (t-r)dr-b p (t)X 



r=0 



Figure 1 shows the limiting default rate F(t) for different values of the con- 
tagion sensitivity /3 . The default rate increases with p c . Fi gure 2 shows 
the limiting default rate F{t) for different values of the parameter a, which 
specifies the reversion speed of the intensity. The default rate is relatively 
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Fig. 1. Comparison of limiting default rate F(t) for different values of the contagion 
sensitivity /3 C . The parameter case is o — 0.9, a — 4, A = 0.5 and Ao = 0.5. 



insensitive to changes in a for shorter horizons; for longer horizons it de- 
creases with a. The limiting default rate is more sensitive to variation in 
the reversion level A, as indicated in Figure 3. Variations in the diffusive 
volatility a of the intensity have little effect on F(t). 

We finally note that the structure of the unperturbed (i.e., f3 c = (3 s = 0) 
dynamics of the intensity (2.1) was crucial in singling out the equation 
(4.3) as the proper macroscopic effect of the contagion (see the proof of 
Lemma 8.2). The calculations in fact hinge upon the explicit formulae for 
affine jump diffusions developed in [10]. In a more general setting we would 
need a more abstract framework (see [16]). 

5. Identification of the limit. We want to use the martingale problem 
(see Chapter 4 of [11]) to show that \i s converge to a limiting process. For 
every / G C°°{V) and /i G E, define 

{f,t*)E = / /(p)/i(dp). 
Jpev 

Let S be the collection of elements $ in B(0 l (V)) of the form 

(5-1) = tp({fl,fj)E, {h,l*)E, (fM,n)E) 
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Fig. 2. Comparison of limiting default rate F(t) for different values of the reversion 
speed a. The parameter case is a — 0.9, j3 c = 2, A = 0.5 and Ao = 0.5. 



for some M £ N, some (p £ C°°(M. M ) and some {fm}rn=i- Then S separates 
£P(V) (see Chapter 3.4 of [11]). It thus suffices to show convergence of the 
martingale problem for functions of the form (5.1). 

Let's fix / S C°°(V) and understand exactly what happens to (f,fi N )E 
when one of the firms defaults. Suppose that the nth firm defaults at time 
t and that none of the other firms default at time t (defaults occur simulta- 
neously with probability zero). Then 

N,n' \N,n' . ^N,n \ , JV.n' 



1 S'^fs' AT V / 



l<n'<N 



i_)i?=— /(P ,\ t )M t + jy/(P A J 



Note, furthermore, that the default at time t means that f* =0 \^ ,n ds = c n , 
so M^ ,n = 0. Hence, 

(5-2) (f,rf)B-{f,&)B = Jijk), 
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where 

N f / iqC 



n'=l ^ ^ ' ^ 

foralli>0, TV G N and n G {1, 2, . . . , N}. 

We now identify the limiting martingale problem for fi N . For p = (p, A) 
where p = (a, A, a, fi c , f3 s ) G V and / G C°°(V), define the operators 

(A/)(p) = ^ ^ (P) " "(A - A) ^(p) - A/(p), 
(£ 2 /)(p) = |£(p). 



Define also 



Q(p) = A/3 c 



for p = (p, A) where p = (a, A, a, /3 , /3 s ) G P. The generator £i corresponds 
to the diffusive part of the intensity with killing rate A, and £2 is the macro- 
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scopic effect of contagion on the surviving intensities at any given time. For 
<1? £ S of the form (5.1), define 

M „ 



m=l 



(5.4) 

X {(£lfm,(J>)E + (Q,V)E(C-2fm,v)E}- 

We claim that .4 will be the generator of the limiting martingale problem. 

Lemma 5.1 (Weak convergence). For any $£5 and <r\ < T2 ■ ■ ■ rj = 
s <t <T and {V>j}/ =1 C B(E), we have that 



lim E 



Proof. For p = (p,A) where p = (a, \,a, /3 C , /3 b ) eV, define 
(£•/)(» = ^ 2 A §TI (P) " «(A - A) |£(p), 



(4/)(p) = a 



2 9A 2( P ) - 7X 9A ( P ) ^' 



x 6 . 



Then £ a is the generator of the idiosyncratic part of the intensity and C b x is 
the generator of the systematic risk. 
We start by writing that 

= + f {A^+A^jdr + Mt, 

Jr=0 

where M is a martingale and 

M 



A 



NA 



m=l 



1 N 

^Ei( £a /™)(^ ,n )+^(^/ m )(pf' n )}M^ 



n=l 



M 



m=l 



X {(£ a /m,/if )i? + e7v(^X t /m,^ V )s}, 



V 



n=l 
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</2, $ )E + Jfrjt), (f M ,^)E + Jfc£t)) 



~ <p((fl,tf)E, (/2,/4 )e, (fM^t )e)}M 
Using Lemma 3.4, it is fairly easy to see that for all / € C°°('P), 



lim E 

N-^oo 



en 



f \(C b x J,^) E \dr 

Jr=0 



0. 



To proceed, let's simplify j£ For each / G C°°{V), t > 0, N G N and 
n<E{l,2,...,N}, define 



jut) ^ x; I(pf' n )^ m " /(^ ^ 

m=l 



Then 



< 



k: 



2.3 



N 2 



d 2 f 



d\ 2 



c 



where K2.3 is the constant from Assumption 2.3. 



Define i(p) = f A for p = (p, A) G P. Setting 

M a 



m=l 



M 



N,n 



(5.5) 



n=l 



m=l 



dx 



x{(Q,^) E (£ 2 f m ,^) E -(Lf,^) E } 



N\ 



we have that 



lim E 

N^oo 



r=0 



0. 



Collecting things together, we have that 



lim E 

JV-»-oo 



A?> 1 dr + I A^' z dr 



,7V, 2 



3=1 



which implies the claim. □ 

We, in particular, note the macroscopic effect of the contagion. 
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Remark 5.2. The key step in quantifying the coarse-grained effect of 
contagion was (5.5). Namely, we average the combination of the jump rate 
and the exposure to contagion across the pool. 

6. Tightness. In this subsection we verify that the family {fi N }]\f^ is 
relatively compact (as a De[0, oo)- valued random variable); this of course is 
necessary to ensure that the laws of the fi N, s have at least one limit point. 
The complication of course is the feedback through contagion. We need to 
show that the system is unlikely to "explode" via feedback. Our calculations 
are framed by Theorem 8.6 of Chapter 3 of [11]; we need to show compact 
containment and regularity of the fi N, s. 

In particular, compact containment ensures that there is a compact set 
K, such that /if will belong to K, for all N G N and t€ [0,T] with high 
probability; see Lemma 6.1. Regularity shows, roughly speaking, that jif — 
is bounded in a certain sense by a function of the time interval t — s, that 
goes to zero as the length of the time interval goes to zero; see Lemma 6.3. 
By Theorem 8.6 of Chapter 3 of [11], these two statements imply relative 
compactness of the family {^ n }ngn in De[0, oo); see Lemma 6.4. 

Let's first address compact containment. 

Lemma 6.1. For each 7] > and t>0, there is a compact subset fC of 
E such that 



PROOF. For each L > 0, define K L = [-K 2 . 3 , K 2 . 3 } 3 x [0,K 2 . 3 ] 2 x [0,L]. 
Then K L CC V, and for each t > and N G N, 



Here K2.3 and ^1^3.4 are the constants from Assumption 2.3 and Lemma 3.4. 
Let's next define 



sup F{pi? iK}< V . 



iVeN 

0<t<T 



E[^(V\K L )] = ±Y,F{\?> n >L}< 



-K"l,T,3.4 

L 



n=l 




these are compact subsets of E. We have that 
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< -K"l,T,3.4 < ^l,T,3.4 

- 2^ (r, 4- ^2/,/774T7 - 2^ 



-(L+j)vvr+7-^(L + i)3/2' 



Since 



the result follows. □ 



We next need to understand the regularity of the // s. For each £ > 
and JVgN, we define 

& t N ^*{\N> n ;0< 8 <t,n€{l,2,...,N}}. 

Let's also define q(x,y) *= min{|a; — y|, 1} for all x and y in M. 

To proceed, let's first consider the L N, s. A useful tool will be the following 
integral bound. Fix T > and suppose that / is a square-integrable function 
on [0,T]. Then for any < s < t < T, 



f f(r)dr<Vt^\ [ T P{r)dr 
Jr=s V Jr=0 



< 1 f y/i 



i(t-.)'/*|i + jT 



Lemma 6.2. Define 

def 1 



TTien E[Ejv] < i{l + i^2,T,3.^} (where K2,t,3.4 ^ ^ e constant from Lemma 3.4) 
and 

E[|Lf - Lf < (t - ^V^E^^jV] 

/or allO<s<t<T. 

Proof. The bound E[Htv] < ^{1 + #2,T,3.4} is clear from Lemma 3.4. 
To proceed, let's write 

t N t 

n=l 
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By the martingale problem for L N , we have that = + Jvl t where M 
is a martingale and where 



N±^Jr=0 



71=1 

Thus, for < s < t, we have (keeping in mind that L N is nondecreasing) 

\L?-L»\=L?-L? = A t N ~A? + M t -M s . 
We then can use (6.1) to see that 

N , t 



A ^- A s<^Ej r _ X r' n dr<(t-s) 



The claimed bound follows. □ 

Of course, P{Lf e [0, 1]} = 1 for all t > and N G N, so compact con- 
tainment (i.e., condition (a) of Theorem 7.2 of Chapter 2 of [11]) definitely 
holds. 

Moreover, by Lemma 6.2 we have that for any < i < T, < u < <5, and 
0<v<5At, 

mL? +u ,L?)\^ t N ML?,Ll v ) < E[Lf +u - Lf \*» ] < S^E[3 N \^ t N }. 

Theorem 8.6 of Chapter 3 of [11] thus implies that {L^jjveN is relatively 
compact. 

Lemma 6.3. There is a random variable Hn with sup^f gN E[ffjv] < oo, 
such that for any <t <T, <u< 5, and < v < 5 A t, 

Proof. We start by using (5.2) to see that 



where 



If ,,N\ _/f,,N\ , a1,N .2,N rjl,N t,2,N 

\J>f*t )e = \J,Mo )e + A t +A t + B t +B t , 



n=1 Js-U 

R 1,JV _ 1 /"* 9/ (t N,n\. I \N,n N.n jjim 



lf' n ) 
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where, for simplicity, we have defined 



.i,JV 



n def 1 



d 2 f, 



+ (-«AT,n(Af'" - Ajv,n) " ^ ( P f } 

Thus, for any < s < i < T, 

+ Etf(Bl> N , Bl> N )\^ s N ] +Etf{B*> N , Bf>»)\*»]} 

< 4{E[|4 lJV - A^H^f] + E[|A^ - A 2 ' N \\J? S N ] 

+ E^ 1 '" - Bl> N \ 2 \&?] + n\Bt' N ~ B 2 > N \ 2 \& S N }}. 

We now need to get some bounds. Due to Assumption 2.3, for any < 
s < t < T, we have that 



\jl n (t)\<UK2- 



Of 



dX 



+ 



c 



This implies 



\A 2 ' N -A 2 ' N \<{K 2 . 3 



Of 



dX 



+ 



c 



l^f-^l; 



thus, by Lemma 6.2 we have that 

E[\A} N -A 2 > N \\^]<{t-s) l /^K. 



2.7, 



Of 



dX 



c 



+ H/II \E[E N \^ S N ] 



for all < s < t < T. To bound the increments of A 1,N , define 



By Lemmata 3.4 and 10.1 we have that supjy 6N E[S^] < oo. By (6.1), 
\Al> N -Al> N \<(t-s)^E[Z%\^]. 
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We next turn to the martingale terms. We have that 
E[\Bl' N -Bl' N \ 2 \3? s N ) 

N „t / ^ , , \ 2 



E 



<E 



71=1 Jr ~ S V 



dr 



^ s 



<(t-s)^E[Ef\^], 
E[\B 2 ' N -B 2 > N \ 2 \JP S N ) 



e 2 N E 



\ 2 




dr 





<e 2 N E 



where 



<s%(t- s f*ElEtf\S,"], 



N r T 



We have that 



f2l (3) 

sup EfHjy ] < oo and sup E[E N ') < oo. 

JVGN JVeN 



Collecting things together, we get that for any < t < T, < u < 5, and 
0<v<5At, 

nq\(f,^ + u)EAf,^)E)\^ t %H(f,^)EAf^lv)E) 
<E[ q 2 ((f,^ +u ) E ,(f,^)EW t N ] 

df 



< 45 1/4 E 



3A 



■11/11 }E N + E$+e 2 N zP 



' □ 



We can now prove the desired relative compactness. 

Lemma 6.4. The sequence {/j, n }n € ^ is relatively compact in De[0, oo). 
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Proof. Given Lemmas 6.1 and 6.3, the statement follows by Theo- 
rem 8.6 of Chapter 3 of [11]. □ 

7. Uniqueness. We next verify that the solution of the resulting mar- 
tingale problem is unique. We will use a duality argument (cf. Chapter 4.4 
of [11]). In particular, here duality means that existence of a solution to a 
dual problem ensures uniqueness to the original problem. 

Lemma 7.1 (Uniqueness). There is at most one solution of the martin- 
gale problem for A of (5.4) with initial condition ir x A D . 

Proof. We will use the duality arguments of Chapter 4.4 of [11]. Define 
E* = |Jm=i C°°(P m ). Let's begin by defining a flow on E* as follows. Fix 
/ G E*. Then / G C°°{P M ) for some M G N. Fix next (pi, p 2 , • • • , Pm) G V m 
where p m = (p m , A m ) and p m = (a m , A m , a m , Pm, P m ) for m G {1,2,. . . ,M}. 
Define 



(r t /)( Pl ,p 2 ,...,p M )= f E 



fiPt »Pt '•••'Pi ) ex P 



m=l- /s =° 



where p t ' = (p, A 4 ' ) and 

A t *' m = A m - a m /" - A m ) + cr m /" V^PdWT 

Js=0 Js=0 

for all m G {1, 2, . . . , M}. We also define 

<9 / 

(H m f)(pi, p 2 , • ■ • , Pm, Pm+i) = M/3 M+1 A A /+i -~t— (pi, P 2, • • ■ , Pm) 

for m G {1, 2, . . . ,M} and p M +i = (pm+1, Am+i) G V where p M +i = ("m+i, 
Xm+x,o-m+i, Pm+v Pm+i)- Suppose that f £ E* and that in fact / G C°°(P M ) 
for some M G N. Let e be an exponential(l) random variable. Set Ft == Ttf 
for t < t. Select m G {1,2, ...,M} according to a uniform distribution on 

{1,2,..., M} and set F e = f H m {T c f). Restart the system. 

Let's now connect F to fi. Fix / G E* and fi G Then / G C°°{V M ) for 
some M G N, and we define 

(7.1) 4>(n,f) d = / /(pi,p 2 ,...,PA/)/u((ipi)/i((ip2)---^((ipM)- 

J(pi,p 2 ,-,PM)e^ M 

If we fix 1 = mi < m 2 < m 3 < • • • < m L+ i = M + 1 and {//}f =1 C C°°(V) 
and assume that 

/( P i,P2,...,Pm) = 1 II /*(Pm)| 

1<1<L mi<m<mi+i — l 
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for all (pi, p 2 , • • • , Pm) G V M , then 



L 

m !+ i-m ; -l 
E 

1=1 



By Stone-Weierstrass, we can thus approximate <!> in S by linear combina- 
tions of functions of the form 0(-, /) of (7.1) for some /'s in E. 

To proceed, let's fix f £ E and apply A to the function // i— > f) given 
by (7.1). It is fairly easy to see that if {fit}t>o satisfies the martingale prob- 
lem for A, then for each / £ E, 

<p(pZ,f)= f hi{KJ)ds + Mf\ 

Js=0 

where is a martingale and where, if / E C°°(V M ), 

M . 

hi{nJ) = Y, / {(A,m/)(p) + <Q J ^( J C 2 ,m/)(p)} 

m=1 ip=(pl, P 2,...,PM)G-P A/ 

x fi(dpi)fj,(dp 2 ) ■ • • ju(dpjvf), 

where £i, m and £2,™ denote, respectively, the actions of £1 and £2 defined 
by (5.3) on the mth coordinate of /. On the other hand, we also have that 
for n e E, 

ft 



[ h 2 {^F s )ds + M?\ 

Js=0 



where A4^ is a martingale and 



h 2 (fi,f)=Y] / (£i )m /)(p)/i(dp 2 )---/i(dpM^ 
m=1 Jp=(pi,P2,...,pM)GP M 



1 M 



m=l 

Note that 

M 



M 

m=l 



M f df 

m=1 Jp=(pi,P2,...,PM,PA/+i)GP M+1 aA « 
X] / (Q,V)E(£2,mf)(p)K d P2)---V(dpM)- 

— 1 Jp=(6i ,62,. ...dm )eV M 
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Collecting things together, we have that 

hx(fJ>, f) = h 2 (n, f) + <p(n, f) 
and this implies uniqueness. □ 

8. Proof of main theorem. We now have our first convergence result. Let 
Qn be the P-law of , that is, 

® N (A) d ^F{» N eA} 

for all A G @(D E [0,oo)). Thus, Q N G &>(D E [0,oc)) for all N G N. For u G 
D E [0,oo), define X t (u) = f w(t) for all f > 0. 

Proposition 8.1. We have that Qn converges fin the topology of 
&(De[0, oo))/ to the solution Q of the martingale problem generated by A of 
(5.4) and such that QXq 1 = 5 nX Ao ■ I 71 other words, Q{Xo = ir x A D } = 1 and 
for all $ G 5 and 0<r 1 <r 2 <---<r J = s<t<T and {ipj} J j=1 C B(E), 
we have that 



lim E^ 

iV-»-oo 



|$(X t ) - $(X S ) - ^' (^$)(X r )dr| JJVi(Xr 



0. 



where E^ zs i/ie expectation operator defined by 



Proof. The result follows from Lemmata 5.1, 6.4 and 7.1. Of course, 
we also have that for any $ G S, 

E®[<t>(X )} = lim E«[$(/^)] = $(ttx A ), 

TV— >oo 

which implies the claimed initial condition. □ 
We next want to identify Q. 

Lemma 8.2. We have that Q = 5p i , where ft is given by (4-5). 

Proof. Recall (4.4) and the operators C\,C 2 from (5.3) and the defini- 
tion of Q in (4.3). For any / G C°°(V), 



(f,fJ>t)E 
Thus, 



E 



E 



p=(p,a)£P 



7(pA*(p))exp 



(A/)(p,A*(p))exp 



X*(p)ds 



s=0 



n(dp)A (dX). 



K(p)ds 



n(dp)A (d\) 
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+ f E 

Jp=( P ,x)eP 



(£ 2 /)(p,A*(p))Q(i)exp 



=o 



x 7r(dp)A (dA) 

= {£lf,flt)E + Q(t)(£ 2 f,p4)E- 
To proceed, define 



G(t) 



dcf 



f3 c E 



'p=(p,A)6P 

p=(a,A,o-, / 9 c ,/3 s ) 



cxp 



A*(p) ds 



s=0 



n(dp)A (dX). 



On the one hand, we have that 



G(t) 



(3 C E 



ip=jp,X)eV 



ir(dp)A (dX) 



p u \fi t (dp) = -(Q,fi t ) E . 

>p=_(p,X)er 

p={a,\,a,(3 c ,fi s ) 

We want to show that 

(8.1) G{t) = -Q(t). 

Indeed, fix p = (p, A) G V where p = (a, A, a, (3 C , (3 s ). Define 



M s = exp 



b p (t-s)X*(p)- b p (t-r){Q(r) + aX}dr- X* r (p)dr 



r=0 



for < s < t. Using the calculations of [10], 

dM s = dM s + {b p (t - s)A*(p) - W(t - s){-a(X* s (p) - A) + Q(s)} 



+ \a 2 (b p {t - s)) 2 X* s {p) + P(t - s)(Q(s) + aX) - X* s (p)}M s ds 



dM s , 



where Ai is a martingale [we use here the ODE (4.2)]. Noting that 



Mo = exp 



b p (t)X- / b p (t-r){Q(r) + aX}dr 



r=0 



M t = exp 



A*(p) dr 



r=0 



we have that 
G(t) = 



(3° exp 



> P ={p,\)eV 

P =(a,X,<r,l3 c ,P s ) 



-b P (t)X 
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b p (t-r){Q{r) + a\}dr 



'r=0 

x 7r(dp)A Q (dA). 
Differentiating this, we get that 



G(t) 



8 



c 



'p=(p,A)eP 



b p {t)X + / b p (t-r){Q(r) + aX}dr 



r=0 



x exp 



b p (t)X- [ b p (t-r){Q(r) + aX} 

Jr=0 



dr 



x ir(dp)A (dX) 



-Q(t), 



where we have used the defining equation (4.3) for Q. Thus, (8.1) holds, so 
we have that 

d -{f,fit)E = {Cif,fk) E + {Q,P-t)E{C 2 f,fit)E- 



dt 



Thus, 



= S(jio)+ / (A$)(p s )ds, 

Js=0 



and, hence, 5p satisfies the martingale problem generated by A. Of course, 
we also have that fio = it x A q . By uniqueness, the claim follows. □ 

We now can finish the proof of our main result. 

Proof of Theorem 4.2. Since weak convergence to a constant im- 
plies convergence in probability, we have (4.6). Using the fact that the map 
ip i *P i — y 1 is in C(P), L N is a continuous transformation of fj, into Dr[0, oo). 
From (4.7) we have that 

lim F{d R {L N ,F)>5} = 

N^oo 

for each 6 > 0. To finish the proof, we need to replace the Skorohod norm 

def 

d-M. by the supremum norm. From (4.8) we have that Kt = sup 0<t<T F(t) is 
finite for each T > 0. To get the claimed convergence, we adopt the notation 
of Chapter 3.5 of [11]. For any nondecreasing and differentiable map g of 
[0,T] into itself and any t€ [0,T], we have that 

\Lf - F(t)\ < \Lf - F(g(t))\ + \F(g(t)) - F(t)\ 

< sup \L? -F(g{t))\+K T \g(t)-t\ 

0<t<T 
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< sup \L? -F{g(t))\+K T T sup \g{t)-l\ 

0<t<T 0<t<T 

< sup \L?-F(g(t))\ 

0<t<T 



+ i^rTmax< exp sup |logg(t)| 

k><t<T 

sup \logg(t)\ -if. 

0<<<T 



Varying g, we get that 



sup \L?-F(t)\ 

0<t<T 

< d R (L N ,F) + K T Tm a x{\exp[ML N ,F)} - 1|, 

-1|}. 

The claim now follows; note that and L w both take values in [0, 1]. □ 

9. Conclusion and extensions. We have developed a point process model 
of correlated default timing in a portfolio of firms, and have analyzed typical 
default profiles in the limit as the size of the pool grows. Our empirically mo- 
tivated model captures two important sources of default clustering, namely, 
the exposure of firms to a systematic risk process, and contagion. We have 
proved a law of large numbers for the default rate in the pool. 

There are several potential extensions of our work. For example, the de- 
fault intensity dynamics (2.1) can be generalized to include a dependence 
on the systematic risk process of the magnitude of the jump at a default. 
Then, the impact of a default on the surviving firms depends on the state of 
the systematic risk: intuitively, if the economy is weak, firms are fragile and 
more susceptible to contagion. This generalization of the intensity dynamics 
is empirically plausible, and can be treated with arguments similar to the 
ones we currently use. 



10. Proofs of Lemmas 3.1, 3.2, 3.4 and 4.1. In this section we prove 
Lemmas 3.1, 3.2, 3.4 and 4.1. For presentation purposes, we first collect in 
Lemma 10.1 some a-priori bounds that will be useful in the proof of these 
lemmas. Then, in Section 10.2 we proceed with the proof of Lemmas 3.1, 3.2 
and 3.4. We mention here that the square-root singularity unavoidably com- 
plicates the analysis. The theory behind CIR-like processes is a bit delicate 
due to the square root singularity in the diffusion, so we need to develop some 
new modifications to existing results (cf. [15, 17, 18]). Last, in Section 10.3 
we prove Lemma 4.1. 
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10.1. Effect of systematic risk. Our first step is to get some usable bounds 
on the systematic risk X . We need these bounds since, as we mentioned in 
Section 3, the XfdXt term contains the term XtXtdt, implying that the dy- 
namics of the M 2 -valued process (A, X) contain a superlinear drift. Note that 
the systematic risk process X of course has an explicit form: 

X t = e-^x + f e -T(*" s ) dV s , t>0. 

Js=0 

Fix p= (a, A, a, /3 , {3 s ) £ V, A in E + , and £ as required in the beginning 
of Section 3. Define 

r t d ^ f at + /3 5 7 f X s ds, 

Js=0 

Z t d = A + aX [ e Vs ds + f [ e r ° d£ s 

Js=0 Js=0 

= X + a~X f e r ° ds + t3 C {e Vt Zt - I e r ^ s (a + flX s ) ds) 

Js=0 I Js=0 J 

= X Q + [ e r °{a~X - f3 c Uu + flX s )} ds + p c £ t e Ft 

Js=0 

for all t > 0. The alternate representations of Z will allow us bounds which 
are independent of £. 

Our first result is a bounds on X, T and Z which explicitly depend on 
various coefficients. The importance of the bound on the moments of Zt is 
that they do not depend on £. 

Lemma 10.1. For each p> 1 and t>0, 

E [^]^<|x | + ^(M) 1/(2P) , 

E[exp[pr t ]] < exp[\p\{at + |/3 c x |} + \{pP C ?t], 
E[z 2 P] i/(2 P ) < Aq + | /3 C| E[e -2 P r t] i/(2 P ) 

/ft \ V(4p) 

+ t i-i/(2p)M E[e -4 P r a]ds \ 

x | a At 1 /(4 P ) +a | /3 C| t l/(4 P ) 

/ ft \V(4p)- 

+ w c f3 s 1 \[j s _nxt p ]ds 
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Proof. We first bound X. For every p > 1 and t > 



E[x 2 P] l/( 2p )^ |xoe _ 74| + | E 



2p- 



1/(2P) 



|x e-^| + 



-2 7 (t-s) rfs 



s=0 



(2p)iy/(2p) 



<\ x \ + J-(M>i 



V(2p) 



Xo + 



1 /(2p)!\ 1/(2p) 



2 V / 7V p! 



Next note that 



T t = at-(3 s j I x e-^ s ds 
's=0 



/3 5 7 /" ( /" e-^ s - r) dV r \ 



s=0 Ur=0 



ds 



at - /3 5 x {l - e" 7 '} - /3 5 7 



r=0 \J s=r 



-y(s-r) ds I dVr 



at - p s x {l - e" 7 '} - (3 s f {1 - e-^- r) } dV r . 

Jr=0 



Thus, for any pGM 

E[exp[^r t ]] = exp 



p{at + p c x (l-e-^)} 



+ 



f {l- e -T(*-r)}2 dr 
Jr=0 



< exp 

We can finally bound Z. We have that 
E [^P]i/(2p) <A +E 



|p|{at+|/3 c x |} + ^(p/3°yt 



1 



C\2j 



t \ 2p 

e" rs {aA - /? c 6(a + P S "fX s )} ds 



s=0 



l/(2p) 



+ |/3 c |E[e-^ rt ] 1 /^. 



We also have that 
r-t 



E 



\ 2p 



•{a\-/3 c t s (a + p s 7 X s )}ds) 



s=0 



V(2P) 



30 



K. GIESECKE, K. SPILIOPOULOS AND R. B. SOWERS 



<E 



<E 



t \ v 

e 2Ts ds 

t 



z 2Ts ds ) 



s=Q 



\ 2p 



{a\- + p s -yX.)} 2 to ' 

s=0 



V(40 



V(2p) 



x E 



{aA-/3 c e s (a + /3 5 7^ s )}^ 

s=0 / 



l/(4p) 



< t 1 ~ 1 /(2 P ) E 



-4 P r a 



r/.s 



l/(4p) 



x E 



{aA-/3%(a + /3 5 7 X s )} 2p ^ 



s=0 



l/(4p) 



< jl-l/(2p) 



E[e- 4prs l(is 



\ V(%0 



V 

x jaAt 1 /^ + al^lt 1 /^) + |^ 7 |E 
Combine things together to get the bound on Z. □ 



Xt p ds 



s=0 



i/(4») 



10.2. Proofs of Lemmas 3.1, 3.2 and 3.4- Let's next understand the reg- 
ularity of various CIR-like processes which we use. Before proceeding with 
the proofs, we define a function ip^ix) that will be essential for the proofs. It 
is introduced in order to deal with the square-root singularity. In particular, 
let 



def 2 



In 77 1 



y 1 



-X[t,,t,i/2](z) dz \ dy and g v (x) d = \x\ - ip v (x) 

y=0 I J z=0 z 



for all xeM. 

Let us then study some important properties of ipr)(x) that will be repeat- 
edly used in the proofs. First, we note that tp v is even, so g v is also even. 
Taking derivatives, we have that 

2 



ip v (x) 



In 77 1 



x 1 

-X M /z](z)dz and i) v {x) = 

z=0 * 111 '/ A 



21 ( \ 

r~XM 1/2 ]W 



for all x > 0. Since tj^ = — ip v < 0, g v is nonincreasing. For x > y/rj, 

, , In77 1//2 — In 77 

8 - W = 1 - 2 Ml/,) - = - 

so in fact 5^ is nonnegative on (0,oo) and it vanishes on [y/rj, 00). Thus, g v 
is nondecreasing and reaches its maximum at y/rj. Since 5r)(0) = 0, we in fact 
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have that 

<g v (x) <g v (Vv) 

for all x > 0. Since g v is nonincreasing on (0, oo) and g v (x) = 1 for x G (0, rj), 
we have that g v (x) < 1 for all x G (0, y^r/), so (^(y^?) < -y^- Since ^ is even, 
we in fact must have that 1^(^)1 < \ft) for ah x G K. Hence, 

M < 1pr,{x) + a/t? 

for all x G R. We finally note that 

1 1 



■■2 1 2 



\X\ 7] 



for all 

Now we have all the necessary tools to proceed with the proof of the 
lemmas. 

PROOF of Lemma 3.1. For each N G N, define 

def [tN\ 



g N (t) 

for all t G [0,T]. For each N G N, define 



N 



Y t N = a 



s=0 



(Y* + Z s ) VOdW: + ? [ ((Y e N N{s) + Z s ) V 0) dV s . 

J s=0 



We will show that (Zt + Y t N )e rt converges to a solution of (3.1) (as N oo). 

As a first step, let's bound some moments. Fix p > 1. For < s < 
t<T, [17], Exercise 3.25, gives us that 

E < vN 



{{Y» N[r) + Z r )V0)dV r 

s 

<{p{2p-l)f{t-sY~ l f E[((Y e N N(r) + Z r )V0) 2p dr] 

J r=s 

<{p{2p-l)f{t-sf~ l f E[\Y e N N{r) + Z r \ 2p ]dr 

J r=s 

<2 2p ~ 1 (p{2p-l)) p (t-s) p ~ 1 

' E[\Y g N N{r) \ 2p ]dr + f E[\Z r \ 2p ]dr\. 

» J r=s J 



Similarly, 



E 



f e rr/2 ^J(Y e N N{r) + Z r )V0dW; 



2p 
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<{p(2p-l)f{t-sY- 1 f E[e^\(Y e N N{r) + Z r )V Of] dr 

J r=s 

< \(p(2p - l)f(t - sf- '{f^n^] + n\Y e N N(r) + Z r \ 2p ] dr 

< l -(p(2p-i)Y(t- s y- 1 

x{T E[e^]dr + 2^ f E[\Y^ [r) \ 2p ] dr 

yJr=s J r=s 

+ 2 2 P -i f* E [\ Zr \ 2p ]dr\. 

J r=s J 

We can bound the effect of Z by Lemma 10.1. Collecting things together, 
and using the fact that f?jv(i) < t, we have that there is a K A > such that 

n\Y g N N(s) \ 2p ]dr 
<K A + K A f E[\Y e N N(s) \ 2p ]dr 

J s=0 

for all N £ N and t G [0,T], which in turn implies that 

(10.1) sup E[\Y e N N(t) \ 2p ]<K A e K - T 

0<t<T y ' 

for < t < T. This in turn implies that there is a Kg > such that 

(10.2) E[\Y t N - Y g N N{t) \ 2p ] <K B \t- QN (tW < K B ± 
for all < t < T. 

We next want to show that Y converges in L 1 . Fix N and N' in N and 
define 

N,N'dei v N V N' 
u t — x t ~ x t ■ 

Fix also 7] > 0. We have that 

I N,N'\^ i , N,N\ , i — 2 »1,N,N' . a 2 a2,N.N' , . A r- 

where M is a martingale and 
1 rt 



A t N ' N ' = ~ 2 j^U^'V'i^Y^ + Z s ) V - ^Y g l {s) + Z.) V 0} 2 ds 
^ ^(^')e r lC( S )-C( S )l^ 

's=0 



1 

< - 

~ 2 
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ri 
's=0 



< 



< 



f J. f r , 1 nWV viV i , 1 i v 



i /■* r r i i ^ 2 



^Fi^ e - + i 1+ -J Y -- Y i:a + -J Y > -Cmi> h» 



r-t 

< 

4 In 77 



A 2 t N ' N ' = \ [ s=o M»?' N '){((Y e N N{s) + Z s ) V 0) - {{Yp {s) + Z.) V 0)}' , 
/ 



3 * 



< - 
~ 2 



s=0 



3 

< 



_/ /r.l/2 , I|vJV_viV |2 l,vJV' _ y A" \2\, 



2 In ?/ 

In the bound on A 1,N ' N ' , we have used Young's inequality, and in the bound 
on 

A 2,N,N> 

we have used the fact that the support of ^ is contained in 
[Ojy^]. Collecting things together, we have that there is a K > such that 



In rj 1 \ Nrj 2 N'r/ 2 J ' 
.2,JV,JV'i < K j 1/2 



1 1 J - lnrr 1 ! iVr? iV'r? 



for all t€ [0,T]. Thus, 



N,N'-+oo lnr/ -1 In?? -1 



for all i G [0, T]. Letting rj \ 0, we indeed get that limjv,jv'->-oo^[| l/ i V ' JV |] = 0- 
We thus have that 

hin" E[\Y t N -Y t N '\) = 0. 

N,N'-i-oo 

For any p > 1, we also have by interpolation and (10.1) and (10.2) that 



hm E[\Y t N - Yf n < hm ^E[\Y t N - Y t N '\m\Y t N ~ Y t N '\^\ 
= 0. 
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Thus, there is a solution Y of the integral equation 

Y t = a [ e r ^ 2 v / (Y s + Z s )V0dW: + (3 [ ((Y s + Z s ) V 0) dV s 

Js=0 Js=0 

such that sup fe [ T i E[|Y t | p ] < oo for all T > and p > 1. Setting Y t = f Z t + Y t , 
we have that Yt € flp>i L p and that 

Y t = Z t + a [ e r ^WY s V0dW: + p [ (Y s V 0) dV s . 

Js=0 Js=0 

We claim that Y is nonnegative. For each rj > we have that 
^(Y t ) X R-(Y t )=ij ri (X )xR-(Xo) + 4- I M Y s)XR-(Y s )e r °/ 2 (Y s V0)ds 



2 

p 2 r t 



s=0 



+ ^ / M Y s)XR-(Y s )(Y s \/0yds + Mt, 

1 Js=0 

where Ai is a martingale. Taking expectations and then letting 77 \, 0, we 
have that E[Yf] = 0. We finally set \ t d = e~ Vt Y t . The claim follows. □ 

Proof of Lemma 3.2. Let A and A' be two solutions of (3.1). Define 

Y t = Aje r * — Z t and Y( = f A^e Fi — Z t . Since A and A' are assumed to be 
nonnegative, Y and Y' satisfy 

Y t = a [ e v °' 2 ^Y s + Z s dW* s +f[ (Y s + Z s ) dV s , 

Js=0 Js=0 

Yl = a f e T */ 2 y/Y> + Z s dW:+/3 s f (Y^ + Z s )dV s . 

Js=0 Js=0 

Set v t = f Y t — Y(. For each rj > 0, 

where M is a martingale and where 

A l = \f Mv s )e r °WY s + Z S - y/Yl + Z s } 2 ds < — ^ /* e r *ds, 

2 Js=0 inr ? Js=0 

Collecting things together, we have that 

Let r\ \t to get that Y = Y' . The claim follows. □ 
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Let's next prove the needed macroscopic bound on the A ' n 's. 

PROOF of Lemma 3.4. For each N G N and n £ {1, 2, . . . ,N}, define 
d f /"* 

#t ' = Aat^o + a N , n \ Ntn / e s ds + / e s dL s 

Js=0 Js=0 



and let Y N,n satisfy the equation 

f-t 



Y t =a N 



n f e^ n '*y/Y a N > n + Z?' n dW? 

Js=0 

+ s N ^ n f (Y s N > n + Z?> n )dV s ; 

Js=0 



then Af ' n = e r i (Y t N ' n + Zf' n ). We calculate that 
|Af 'T < \{e~ 2 ^ + \Y t N > n + zf' n | 2p } 

< l{ e -2prf- + 2 2 P -l^ Y N,nfp + \Z t N > n \ 2 P)}. 

From Lemma 10.1, we have that 

N N 
-P ^Y,ne 2pT ^] and sup lj>[|Z^] 

0<t<T ^ TZ-y 0<t<T W ^ x 

NeN NeN 

are both finite. 

For each N 6 N and n £ {1,2, ... , N}, we compute that 

n\Y t N ' n \ 2p ] =^-i){^ n £ o E[|rf.n^- 2 e r ^|r^" + ^|]^ 

+ 4(/4.J 2 T E[|y^|^- 2 |y^ + zf'«| 2 ]d s 

Js=0 

To bound the integrals, we have that 

\Y s N > n \ 2 P- 2 e r °' n \Y s N ' n + Z?> n \ 

< jL e 2prf' n , P~ l |yjy,n.ap , J_|yiV,n , z N, n] 2p 

~ 2p P 2p 

< J-paprf'" , P~ 1 iyJV,ni2p , 2 ^ 1 riyAr.n|2p , | Z 7V,n,2 n 
|y7V,n|2p-2|yAf,n _|_ ^JV,n|2 
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1 



< 



< 



P ^|yJV,ni2p _|_ t\y N ' n + Z N,n \ 2p 



P 

P ~ 1 |yJV,w|2p 



P 

+ - {|F/^ + |Zf> n | 2p }- 



Combining things together, we have that there is a -fT > such that 

E[|y^|^<^{^ in + 4(^,J 2 } 

rt 



E[\Y A N ' n \ 2p ]ds+ [ E[e 2pr "' n ]ds 

Js=0 

E[|Zf< n | 2p ] ds 



t 

s=0 



s=0 

for all N G N and n G {1,2, .. . , N}. Using Assumption 2.3 and averaging 
over n, we get the claimed result. □ 

10.3. Proof of Lemma 4-1- Define a homeomorphism $> of C[0,oo) as 

r-t 



*(?)(*) 



clef 



/3 



c 



>p=( P ,\)eV 

P ={a,\,(7,l3 c ,/3 s ) 



b P (t)X 



b p (t — r){q(r) + aX} dr 



r=0 



x exp 



-b p (t)X- / fe p (t-r){g(r) +aX}dr 

Jr=0 

x 7r(dp)A ((iA) 

for all q G C[0,oo) and t > 0. Note that since 6, g and A are all nonnegative, 



< exp 



-b p (t)X 



b p (t — r){q(r) + aX} dr 



r=0 



< 1. 



We can then set up a recursion; we want to solve Q = & (Q). Note that there 
is a K > such that 



\$(q)(t)\<K q(r)dr 
Js=0 

for all nonnegative q G C[0,oo). 

For any q\ and 172 in C[0,oo), we have that 



$(?i)(t) - $(<&)(*) =r t (gi 5?2 ) + r t b (gi,g 2 ), 



where 
r?(gi,ga 

def 



/3 c 6P(t- S ){ 91 ( S ) - 92 ( S )} 



=0 1.7(9=0 Jp=(p,A)€P 

p=(a,A,cr,/3C j /3S) 
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x exp 



-b p (t)X 



r=0 



b P (t ~ r)[{q 2 {r) + %x(r) - q 2 (r))} + aA] dr 



x 7r(dp)A (dA) d0 > ds 



r?(gi,«2) 

def 



s=0 



/3 c <{6 p (i)A+ / feP(t-r) 



p=(P,A)GP 



r=0 



p=(a,A,o-,/3 c ,/3 s ) 

x [{g 2 (r) + %i(r) - 92 W)} + aA] dr 
x{6P(t- s )( gi ( s )-g 2 ( s ))} 



x exp 



-6 p (t)A 



& p (* " r)[{<fe(r) + 0(gi(r) - g 2 (r))} + aA] dr 

r=0 

x 7r(dp)A (dA)d6'| ds. 

Standard techniques from Picard iterations give us the result. 
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